Magnetic field driven instability in planar NJL model in real-time formalism 
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It is known that the symmetric (massless) state of the Nambu-Jona-Lasinio model in 2+1 di- 
mensions in a magnetic field B is not the ground state of the system at zero temperature due to 
the presence of a negative, linear in \a + in\, term in the effective potential for the composite fields 
a ~ {jjtb and tt ~ ipi~/ ip, while the quadratic term is always positive (a tachyon is absent). We 
find that finite temperature is a necessary ingredient for the tachyonic instability of the symmetric 
state to occur. Utilizing the Schwinger-Keldysh real-time formalism we calculate the dispersion 
relations for the fluctuation modes of the composite fields a and tt. We demonstrate the presence of 
the tachyonic instability of the symmetric state for coupling constant that exceeds a certain critical 
^ ' value which vanishes as temperature tends to zero in accordance with the phenomenon of magnetic 

, catalysis. 

PACS numbers: 11.30.Qc, 11.30.Rd, ll.10.Wx 
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^ 1 i. introduction 

i 1 For many years relativistic quantum field models in (2+1) dimensions have attracted a significant interest both 

, due to their sophisticated dynamics and the fact that they describe long wavelength excitations in several planar 
'j^l condensed matter systems [l[ , among them graphene @ , the d-wave state of high T c superconductors Q , topological 
O I- insulators Q and optical lattices [H . Recently there has been a surge of activity in this area connected with the 
experimental discovery of graphene (| whose quasiparticle excitations are described by the massless Dirac equation 
in (2+1) dimensions that leads to many unusual electronic properties of this material and opens new perspectives for 
electronic devices (see, review papers [7[). Lattice effects necessarily produce local interactions for quasiparticles in 
y—i [ graphene [|| and, thus, one naturally comes at the gauged Nambu-Jona-Lasinio (NJL) model in 2+1 dimensions. 

Historically the Nambu-Jona-Lasinio model Q was the first model in which the mass generation and dynamical 
' symmetry breaking (DSB) were considered in elementary particle physics and quantum field theory. At present 
NJL-type models have a significant practical value, for example, the NJL model provides a successful effective theory 
' of low energy Quantum Chromo dynamics [Iol - [l2T ] . Dynamical symmetry breaking occurs in the NJL model only in 
• , supercritical regime when its coupling constant G exceeds a critical value G c . This is different from the Bardeen- 
Cooper-Schriffcr (BCS) theory where a gap in quasiparticle spectrum is generated for any value of coupling constant. 
The physical reason for zero value of the critical coupling constant is connected with the presence of the Fermi surface 
in the BCS theory. According to the renormalization-group studies (l3j . the renormalization-group scaling takes 
place only in the direction perpendicular to the Fermi surface that lowers effectively the space-time dimension by two 
units to a (l+l)-dimensional theory, where as is well known, symmetry breaking occurs for arbitrary weak attraction 
between fermions. 

Since dynamical symmetry breaking in (3+1) and (2+l)-dimcnsional theories requires strong coupling (g c > 1) it 
makes the quantative study of DSB a difficult problem. Therefore, it is very interesting to consider field-theoretical 
models where DSB takes place in the regime of weak coupling (g c w 0). The DSB in a magnetic field [3, EH (magnetic 
catalysis) gives the corresponding example (see also Refs.[l6| and a short review Ref.[T3|)- The essence of the magnetic 
catalysis phenomenon is that the dynamics of the electrons in a magnetic field, B, corresponds effectively to a theory 
with spatial dimension reduced by two units (note a close similarity with the role of the Fermi surface in the BCS 
theory) if their energy is much less than the Landau gap The zero-energy Landau level has a finite density of 

states and this is a key ingredient of magnetic catalysis which plays, in fact, the role of the Fermi surface. 

The magnetic catalysis is an universal phenomenon and its main features are model independent fbH Il8j . It was 
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studied, besides a (2+l)-dimensional NJL-type model, in the NJL3 + i model [15|, quantum electrodynamics [19(, and 
quantum chromo dynamics [20| . The universality of this phenomenon is confirmed by applying holographic techniques 
which have proven to be a powerful analytic tool in studying the qualitative properties of strongly interacting physical 



systems such as interacting quark gluon plasma, graphene, superconductivity, and superfluidity [21 1 . 

In the theory of superconductivity the normal state of a superconductor is unstable at sufficiently low temperature 
with respect to the transition to a superconducting state. This instability is signaled by a pole in the scattering 
amplitude of the electrons with opposite momenta and is known as the Cooper instability 22 1. This instability is 



resolved in the superconducting state through the formation of a condensate of Cooper pairs 



23| that opens a gap in 
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the electron quasiparticle spectrum. 

The instability of the normal state of a quantum statistical or field system has a precursor in the corresponding 
one particle problem which is known as the fall-into-the-center phenomenon. For example, in the study of dynamical 
chiral symmetry breaking in strongly coupled QED [241 ] . the corresponding one particle problem is formulated as the 
Dirac equation for the electron in the field of the Coulomb center and the precursor of the normal state instability 
in QED corresponds to the supercritical charge problem when the lowest in energy bound state dives into the lower 
continuum. Then an electron-positron pair is spontaneously created from vacuum with the electron shielding the 
supercritical charge and positron emitted to infinity (described by a resonance state) [25l , |2(| . 

It is interesting to see what is a precursor of the magnetic catalysis phenomenon in quantum field theories and what 
are its characteristics. Recently the corresponding study was performed in the case of graphene in Ref.[27j, where the 
Dirac equation for the electron in the field of the Coulomb center in a magnetic field was considered and it was shown 
that, as suggested by the magnetic catalysis phenomenon, indeed any charge in the gapless theory is supercritical. 
However, no resonance state was found that is related to the fact that charged particles cannot propagate freely to 
infinity in a constant magnetic field in two dimensions. Still it was found that the low energy bound state crosses 
the level of filled states that suggests that the normal state of the system in a magnetic field should suffer from a 
tachyonic instability (i.e, an analog of the Cooper instability in the theory of superconductivity should exist). 

In the present paper, we directly address the problem of instability of the symmetric state of quantum field theories 
with attraction between fermions and antifcrmions in a magnetic field in the framework of the NJL2+1 model. The 
model is described in Sec. [TTJ The analysis of the effective potential indicates the necessity of finite temperature for 
the tachyonic instability to be present. In Sec. using the Schwinger-Kcldysh real-time formalism, we calculate the 
dispersion relations for composite fields in the LLL approximation and for sufficiently low temperature find a tachyonic 
instability. The contribution of higher Landau levels to the dispersion relations for composite fields is considered in 
Sec. II VI The main results are summarized in Conclusion. 



II. MODEL AND EFFECTIVE POTENTIAL 



The NJL action in (2+1) dimensions in a magnetic field reads 



fyrfDrf + % [ (^(x)V(x)) 2 + {${x)i^i>{x)f 



(1) 



where = <9 M + ieA^ with the vector potential = (0,Bx,0) that describes a constant magnetic field in the 
Landau gauge. We use four-components spinors corresponding to a reducible representation of the Dirac algebra like 
in Ref. [141]. According to magnetic catalysis, we expect that the symmetric state of model ([T]) is unstable for any 
Go > 0. In order to see this, we calculate the effective potential for the composite fields a ~ 'tpip and it ~ ^i^tp. 
Using the Hubbard-Stratonovich method of auxiliary fields, model (TTJ) can be equivalently rewritten as follows: 



d 3 x 



2G 



(2) 



Assumin g th at a = const and ir = const, the effective potential for composite fields a and 7r was found in the second 
paper in |14| (for more details of calculation see Ref. |28j|). The following propagator for fermions with mass m = (a) 
was used in the derivation: 



G(x,x') 



= i<£(x,x') 



G(x-x'), 



(3) 



where the Schwinger phase [29[ is separated from the translation invariant part G(x — x'). The translation invariant 
part of the propagator can be expanded over the Landau levels (compare with Ref. [1J|) an d in the mixed cj,r 
representation it has the form 



G(w;r) 



2irV 



■ exp 



4/2 ) 



El 



(70W + to) 



7r£j,_i 



2P 



(4) 



where V± = (1 ± isgn(ei?)7i72)/2, E n = yjm 2 + 2\eB\n are the Landau levels energies, 1 = 1/ \AeB\ is the magnetic 
length, functions L"(x) are the generalized Laguerre polynomials, and by definition, L n (x) = L^(x), L^_ 1 (x) = 0. 
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Further, according to [l4|, HH, Gil: the lowest Landau level (LLL) contribution is responsible for zero value of the 
critical coupling constant. Since we are interested in the instability of the normal state of the model in the weak 
coupling regime, it is clear that only the dynamics in the LLL can produce this instability. Eq.Q implies that the 
LLL fermion propagator in momentum space is given by 

2,2 V- 

Glll(u;p) = ie p — T —— . (5) 

The effective potential for composite fields a and tt at zero temperature and zero chemical potential in the model 
under consideration was calculated in Ref. 1141. 



V(p) = - 

7T 



A (sfr \ 2 V2 r ( 1 (pi) 2 , A P 

1 )p --s-C -o;-5- + 1 



2y^ V 9 J 13 V 2 ' 2 J 21 



(6) 



where p = \/ o 2 + it 2 , g = GqA/it, A is the UV cut-off, and £(s, q) is the generalized Ricmann zeta function. For 
p — > 0, at weak coupling g <C ypn we have 

The presence of negative linear term p/(2tt1 2 ) clearly indicates that the true minimum of the effective potential 
corresponds to a state with broken symmetry. However, the second derivative of the effective potential with respect to 
p is always positive, hence the tachyon is absent. This situation is rather unusual and the reason for the existence of 
the linear term was explained in )l4| . Since the gap equation is given by dV(p)/dp = 0, the effective potential can be 
reconstructed up a constant by integrating the gap equation. It suffices to consider only the field a (the dependence 
of the effective potential on 7r can be easily restored using the chiral symmetry). Further, the gap equation in the 
model under consideration equals 

a = G (0\W>\0) ■ (8) 

The point crucial for the existence of the linear term in the effective potential is that the chiral condensate (O^^IO) 
does not vanish as a — ¥ even in the free noninteracting theory. It suffices to keep only the LLL contribution. Then 
we have 

(0|^|0) = -i lim -—3 / dujd 2 p ^ = • ( 9 ) 

Integrating it, we find the linear term in the effective potential (JT]). Since fluctuations are described by quadratic terms 
in the effective action, the linear term present in the effective potential does not influence fluctuations. Consequently, 
a tachyonic instability is absent in the spectrum of the model at zero temperature. 

The situation changes crucially at finite temperature. At T ^ the effective potential was calculated in Refs. 
[3, H^]. For small p/T -c 1, at weak coupling it is given by 

Vt( P )= {^r- „ ,1 G c (T,B)=4nTl 2 . (10) 



Go G C (T,B) J 2 

Clearly, for Go > G C (T, B), we have an instability of the conventional (tachyonic) type. The critical coupling G C (T, B) 
tends to zero as T — > 0, and the symmetry broken ground state occurs at arbitrary small attractive interaction in 
accordance with the phenomenon of magnetic catalysis. The absence of the linear term at finite temperature is 
consistent with the absence of chiral condensate in the free theory at T ^ 0, 

+ 00 _ 2g2 

(0|#|0} T = - lim ^ / d 2 P (7rr(2? :; P 1))2 + (j2 - - Hm ^ tanh £ = . (11) 



This result suggests that in order to find a tachyonic instability, we should study quadratic fluctuations of the model 
at finite temperature. To do this, we will calculate in the next section the dispersion relations for composite fields a 
and 7r at finite temperature in the LLL approximation and analyze them. 
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III. TACHYONIC INSTABILITY FOR COMPOSITE FIELDS IN THE REAL-TIME FORMALISM IN 

THE LLL APPROXIMATION 



The analysis in the previous section shows that in the model under consideration a tachyonic instability can appear 
only at finite temperature. Since instability is an inherently dynamical process and the Matsubara imaginary time 
formalism is mainly used for the study of theories at thcrmodynamical equilibrium, in order to analyze the tachyonic 
instability we will utilize the Schwinger-Kcldysh real-time formalism [3(1 EJ (for a review, see [32[). The action in the 
real-time formalism contains integrals over positive time branch t+ and negative one t— . Then action ([5]) transforms 
into 



S = / d 3 x ip(i^D^ -a - mf)i> - / d 3 x 



2G 



where the time integration proceeds along the closed path time contour [3j 



dt = 



+OC 



dt_ = 



+00 



+00 



dt A 



+00 



dt. 



(12) 



(13) 



Since the fields a± and n± are defined on the positive and negative time parts of the contour, in what follows it is 
convenient to consider their linear combinations 



c c , A 



7T-L ± 7T_ 



7T C . A = 



2 2 

Integrating over fermions in the functional integral, we find the following effective action for the composite fields: 



(14) 



Scir — — 



d 3 x- 



2G 



iTrln [G*- 1 ] , 



(15) 



where G~ l = —{(ij^D^ — a — iir'y 5 )S p (x — y) and S p (x — y) is a contour <5-function, the trace Tr in Eq. (|15j) is taken in 
the functional sense. The equations of motion for composite fields follow from this effective action and are given by 
(the physically sensible case corresponds to oa = ^A = 0) [HI, [33| 



2a c (x) 

Go 
2tt c (x) 



.JTrln [G- 1 ] 



Scta(x) 
.dTrln [G*- 1 ] 



<TA="A=0 



-iTr 



-iTr 



5G- 



6(ta(x) 
5G- 1 



6ita(x) 



Go Stta(x) 
Here G{x, y) is the two-point correlation function defined as 

G(x,y) = Tr(T p (^(x)j(y))p) = (T p (ijj(x)$(y)) 



CTA=TA=0 



(16) 
(17) 

(18) 



p is the thermal density matrix and T p is the time-ordering operator along a complex path p. Since x, y can take 
values on either positive or negative time branches, it is convenient to represent G(x 7 y) as 2 x 2 matrix: 



G(x,y) = 



G + - 



G+- 
G— 



( (T(i>{x)i>{y)) -${y)Tp(x)) 



(19) 



where T and T are the usual time-ordering operator and anti-time-ordering operators, respectively. Note the identity 

G++ + G = G^ + G |- which follows from the identity for the step functions 8(x — y) + 8(y — x) = 1. 

Since we are interested in physical excitations, we will consider time dependent solutions of the above equations 
which deviate weakly from constant values a c (x) = cr{x) + a and 7r c (a:) = 7r(x), where a = const. Then we obtain 



2(g(x) +a) 
Go 



= — zTr 



G 



5G- 



Scta(x) 



(7 A— 7TA— 0, <y c 



+ i d 3 y Tr 



q SG^ 1 Q 5G- 1 



So-a{x) Sa c (y) 



O"A=7l"A=0, cr c 



*(V) , (20) 



27r(a;) 



= — iTr 



G 



5G- 



5cta(x) 



fTA- 7TA— 0, <7 C — G 



+ i d 3 y Tr 



c 6G- 1 G 8G- 1 



Sir A (x) 5ir c (y) 



(7A— 7TA— 0, (7 C — O 



(21) 
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In the single time representation S p (x — y) = t 3 5(x — y), where T3 is the third Pauli matrix, and the variational 
derivatives and functional traces are calculated according to the rules 



SG 1 (x,y) . Sa{x) . . 

= l s aA ( z } °p( x ~v) = l °p\ x ~ vpp\ x - z ) = l °( x - y) d i x - z )i 



(22) 



Tr 



G 



5G- 



d 3 ud 3 v tr 



d 3 ud 3 v tr 



G(u, v) 



t 3 G(u,v)t 3 



8G- l {v,u) 
8<7 a(x) 



Scfa(x) 



= I d 3 ud 3 v tr 



G(u, v 







(23) 



To calculate the right-hand sides of Eas.(|20p and (f2"Tj) , we should determine the fermion Green's function in the real- 
time formalism. As wc discussed in the previous section, for our purposes it suffices to use the LLL approximation. 
The LLL Green's function in the real-time formalism equals 



/ G - * - (t — t') G^ (t — t') 
G LLL (x,x ) - / _A(x,x ) ( G -+ {t _f ) G—(t-t) J ' 



(24) 



where 



^ (X ' X ' )= 2^ eXP 



(x - yi' f 
Al 2 



i$(x, x') 



is the space dependent part of the LLL fermion propagator and 



G++H = i ^X~° - - 2^(7°^ + a)n F {a)5{u 2 - ° 2 ) ■ 
ui z — a A + ic 



(25) 



(26) 



G (w) = -i J^-l —■ ~ 27r (7° w + o)n F {d)5{io 2 - a 2 ) . 



(27) 



G+-(uj) = -2tt ( 7 °w + a) n F {uj)sgn(uj)S(uj 2 - a 2 ) , 



(28) 



G"+(w) = 2ir (j°uj + a) n F (-uj)sgn(Lj)S(uj 2 ~ a 2 ) 



(29) 



are the Fourier transforms of G y (t),i,j = ±. Here n F (a) = (exp(cr/T) + l) 1 is the Fermi-Dirac distribution function. 
Further, it is convenient to perform the unitary Keldysh transformation (3ll . l32l | 



where 



G 



U^GU 



G a 

G r G r 



u 



-( 1 1 , 

V2 V- 1 



G„(w) 
G r (w) 
G c (w) 



- (G++ - G + G+- - G~ 



(G H 



G" 



G^ 



G~ 



G++ - G" 
G++ - G^ 



7 u w + a 



uj 2 — o 2 — iesgnw 
7°oj + (7 

— a 2 + iesgnw 



G ++ + G~ = G+" + G"+ = 27rtanh — ( 7 °u; + a) S(w 2 - a 2 ) 



(30) 

(31) 
(32) 
(33) 



are the advanced, retarded, and correlation functions. For time dependent and spatially homogeneous modes a(t) 
and 7r(<), Eqs. (j2"0|) and (|2"Tj) imply the following equations: 



a(t) =G / dt'W(t-t')a(t'), 



(34) 



= Go / dt'U^(t-t')Tr(t'), 



(35) 
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where 



n CT (i - t') = ^tr [G r (t - t')G c (t' - 1) + G c {t - t')G a {t' - 1)] = ; 



(36) 



W{t - t') = — ^tr [i7 5 G r (t - t r )ij 5 G c {t' - t) + t 7 5 G c (t - t')ij 5 G a (t' - tj\ = tanh ■ ff 



ttZ 2 



2T / 2tt 4ct 2 - ft 2 



(37) 



Hence Eq. ([34| gives <7(i) = 0. Note that the equality IT 7 ^) = is due to the LLL approximation used in this section. 
On the other hand, W(t) ^ in the same approximation. In the next section we obtain expressions for W^it) where 
all Landau levels are taken into account. 
For the Fourier transform 7r(fi), we find 



2nl 2 4a tanh ^ 



Go 



4ct 2 - n 2 



?f(f2) = 



(38) 



that implies 



fl 2 = 4a 



Go a 

a tann — 

2tt/ 2 2T 



(39) 



For a ->• 0, 



n 2 = 4a 2 



Go 



4nTl 2 



(40) 



Obviously, for T less than the critical value 



Tr 



Gq 

I^z 2 



(41) 



we have a tachyon. This result is perfectly consistent with the effective potential at finite temperature (jlpp whose 
symmetric and symmetry broken phases are separated by the curve 



1 



1 



Go 4tt1 2 T 



0, 



(42) 



that leads to the critical temperature (|4"Tj) . 



IV. THE EFFECTIVE ACTION FOR COMPOSITE FIELDS IN THE REAL-TIME FORMALISM 

BEYOND THE LLL APPROXIMATION 

In Sec. Ill, we calculated the correlators 11°' (t — t') and II 7r (< — t') given by Eqs. (f3l>j) and (f3"T|) in the LLL approxi- 
mation. In the present section, we calculate these quantities taking into account the contribution of all Landau levels. 
In addition, we determine the dependence of W and IF 1 ' on spatial coordinates, i.e. calculate II' 7 (t — t' ,x — x') and 
IT 7r (i — t' , x — x') (note that W and IF are translation invariant in spatial coordinates because the Schwinger phases 
cancel out for a closed fermion loop with two vertices). 

In the real time formalism the propagator in a magnetic field and at finite temperature can be written in the form 

<**.»> = I ^-<.°-«°>»giwx- y ,„) ( g:j:> ) . <«> 

where the factor X(x,y) is given by Eq.([25|) and 

D B (r,u,) = tf. + a) (v.L n (^) +V + L n ^ (J)) -i^L^ (^) . (44) 
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Further, 



G~-(oj) = - LJ 2_E2_ l0 - 2™ F (E n )6(uj 2 - £ 2 ), (46) 
G-+(w) = 2™ F (- W )sgnH<5( W 2 - G +- (w) = _ 2mF (w)sgn(u;)5(a; 2 - £ 2 ), (47) 



and = ^ct 2 + 2\eB\n. 

The correlators II cr (a; — x'), IT T (x — x') (x = (i,x)) are defined by the expressions, 

n CT (x - x') = -iTr [G(x, x')G{x', x)] , W{x- x') = -iTr [i7 5 G(x, x')«7 5 G(x', a:)] , (48) 

and the trace Tr includes also the trace (tr) over Dirac indices. Performing the Keldysh transformation the matrix in 
Eq.(g3J} takes the form like in Eq.([30]) with 



G ria {u,n)= a * = {l ,.\ 2 „* > G c ( W ,n)=27rtanh|^ 2 -£ 2 ). (49) 

Taking the Fourier transform of Ii a (t, x) we obtain 

IT(0,k) = -ij^- J d 2 re l]ir tr[G r (r,L})G c (~r,L} + {}) + G c {r,uj)G a (-r,oj + 
°° f dui 

= — i > / — [G r (w, n)G c (uj + fi, m) + G c (oj, n)G a (u) + SI, m)] 

— ' J 27T 

n,m— 

X / ^0^^ r ~ r2/2l2t ^ D n(r^)D m (-r,uj + n)}. (50) 

The space integral equals 

r re 4kM2/2 ' 2 tr[A l (r,w)Ar l (-r,w + 0)] = 2s nm (y)(uj{uj + SI) + a 2 ) - 4\eB\r nm (y), (51) 



2nP 



where 

Sn m (y) = J ^-' /212 (A. (^) ^ (£) +in-! (^) V2/2 

= ^4 e "M £ m m (») £ r n (tf) + ^-™(v)i""i n (»)) - 2/ = k 2 / 2 /2, (52) 

r nm (y) = J ^ ^ r - r2/2l2 Ll-x (^) Ll_ x (^j = (-ir+ n e-y m L^ m (y)L^r(y) • (53) 
[for the evaluation of the integrals s nm (y), r nm (y) see Appendix A in Ref.[3J|]. Therefore, we get 

((w - Q)w + CT 2 )s nm (y) - 2\eB\r nm {y) 



J ~ , n 



n,m— 
2> 



(w - n + «o) 2 - # 2 

((w + Q)uj + (T 2 )s nm (y) - 2\eB\r nm (y)~ 



E 



(w + - iO) 2 - £ 2 

f m ggXg| - £g - » 2 )gnrn(y) - (jg ~ g| ~ ^(^(^CX 2 ~ 2|eB|r nm ( g )) 

E m nl 2 {(E m + E n y-W)((E m -E n ) 2 -W) 



+ (m O n) 



(,54) 



where i m = tanh(-E m /2T). The calculation of the correlator 11" results in the same expression except tr 2 is replaced 
by —a 2 . Thus, wc write 



ir>*(fi,fc)= E 

n , m — 



t m - g 2 - » 2 ) S „ m (y) - {El - Ej - n 2 )(±s nm (y)a 2 - 2\eB\r nm (y)) 

E m 7rP((E m + E n y - a 2 )((E m - E n y - r> 2 ) 



(m -H- n) 



(55) 
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where ± signs correspond to W and IP r , respectively. 

To find the dispersion laws at small and |k| it is convenient to evaluate the sum over the Landau levels. This 
can be done explicitly if temperature is much lower than the value of a magnetic field, T <C y/\eB\. The details of 
calculations are given in Appendix. The dispersion relations are given by the equations, 

-^r+U^(Q,k)=0. (56) 
G 

For a ^ 0, the dispersion relations for a and tt modes at small fH <C 1 and y = k 2 l 2 /2 <C 1 take the form, respectively, 
1 V- 2 IQ- ( P- a(l-tanh(q/2r)) ^ 

J_ V + tanha/2r-l o2 /_ZQ+ tanh(<x/2r) - 1 \ / P+ 1 - tanh(a/2T) \ _ n 
Go \ 3 / 2 l 27r/ 2 a + ^ 3/2 / U 3/2 « 2 ^ 2 * / 

where the quantities V /± ,Q ± ,P ± are given by Eqs. (|A16|) - Eqs. (|A19|) . At the minimum of the effective potential the 
7T mode corresponds to a Nambu-Goldstone boson and a satisfies the gap equation, 

1 V+ tanh £p ~ 1 

-G- + ^m + ^^-°- (59) 

The gap equation written in the form (£(s, v) is the generalized zeta function), 

-^ + 1*^^ + ^(1,1 + ^1!) =0, S = *(±--±-), 6V = ^, (60) 

is in agreement at T = with the one obtained in Ref . [Til . Fixing the intrinsic scale S the gap equation determines 
a as a function of temperature T and magnetic field eB = l/l 2 [We recall that in the used approximation Tl <C 1]. 
The critical line separating symmetric and symmetry broken phases is obtained from Eq. (|60[) when a — > 0: 

^ = 2M-V2CQ), (61) 

and in the weak coupling limit Go <C Gq c it agrees with Eq. (|4"2"|) . The gap equation was analyzed in Ref.[l4j 
at T — where three regions of different behavior of a as a function of a magnetic field were revealed. In the near 
critical region Go — Gq c , where \S\l <C 1, the dependence on the ultraviolet cutoff A disappears and we find al ~ 0.45. 
In other two regions, subcritical (Go < Gq c ) and supercritical (Go > Gq c ), and for \5\l ^> 1, the solution of the gap 
equation behaves 

al ^ ^ < 1, 6 > 0, (62) 
al ~ \S\l > 1, S < 0. (63) 

At finite temperature a nontrivial solution for ct in subcritical region (S > 0) exists for magnetic fields satisfying 

1 



IS < 



C 2 (l/2) + ^ 



2^2 

i.e., for magnetic fields exceeding some critical value. The dispersion laws take the following form: 



(64) 



" 2 = ^ 2 , vl = - 

V2 tanh ^ + (ctZ) 3 C ( f , 1 + ^ 



2T 

for the 7r mode, and 

2 = A/ ff 2 + ^ k2j (66 ) 
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FIG. 2: M%/5 2 dependence on dimensionless magnetic field Z|<5| for temperature T = 10 5 |<5|. Solid (dashed) line corresponds 
to 6 > (8 < 0). 



8 V2tanh^ + (a0 3 C(|,l 



M = ^ '- — (67) 

° ^2c(i,i + ^)-(^) 2 c(ii + ^ x: 



V2al tanh £ + (al) 2 ( ( |, 1 + ) + I |, 1 + M 



- i\ 2 



"2 = 4 ~ 7^ (68) 



c(|,i+^)-^c(f,i 



for the a mode, respectively. At zero temperature Eqs. (j6"5)) - (j6"5j) coincide with those obtained in Ref.pjj]. One can 
check that the quantities v 2 ,v 2 are positive and remain always less than 1 (we set the velocity of light c = 1). Their 
behavior for a chosen value of a temperature (T = 1CP 5 |(5|) is shown in FigQ] where the gap ct(T, /) is determined from 
Eq. (|5Ul) . The behavior of M 2 /S 2 as a function of the magnetic field is shown in FigJ5] All dimensionful quantities in 
Figs lll2l are measured in units of \6\. Asymptotical behavior of the quantities v 2 ,v 2 in subcritical and supercritical 
regions and for \Sl\ ^> 1 is given by the expressions (for simplicity we take T = 0), 

S> °> v 2 ~{ COT' S> °> (69) 



and for the square mass M 2 , 



-W' 5>0 ' CT li-w. 5>0 ' 



cf^ji; 5>0 ' (70) 

* 6^ 2 (l+2W ))• *<°- 
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These asymptotics should be compared with those obtained in Ref . fl4| . 



V. CONCLUSION 

According to the magnetic catalysis phenomenon [l4j |. an arbitrary weak attraction between fermions and their 
antiparticlcs leads to chiral symmetry breaking and gap generation in (3+1) and (2+l)-dimensional theories in a 
magnetic field. Consequently, the normal state of these theories should be unstable in a magnetic field even in the 
weak coupling regime. It is worth noting that the instability of the normal state in a magnetic field is qualitatively 
different for theories in (3+1) and (2+1) dimensions. Since constant magnetic field effectively reduces [HI, HH the 
spacetime dimension by two units for fermions in the infrared region, (3+l)-dimcnsional theories are reduced to 
effective (l+l)-dimcnsional theories, where bound states are easily formed in the weak coupling regime and resonance 
states describing emitted antiparticlcs propagating to infinity are realized in the standard way in the corresponding 
quantum mechanical one particle problems. 

As noted in Introduction, the situation is different in (2+l)-dimcnsional theories. The dimensional reduction in a 
magnetic field means that the corresponding effective theories are (0+l)-dimensional ones. Consequently, no emission 
to infinity is possible. This conclusion is explicitly confirmed by the study of the (2+l)-dimensional Dirac equation 
for the electrons in the field of the Coulomb center in graphene in a magnetic field performed in Ref. [27| where no 
resonance state was found. 

In the present paper, in order to study the normal state instability connected with the magnetic catalysis phe- 
nomenon in a (2+l)-dimensional theory, we considered the weakly coupled NJL2+1 model in a magnetic field at 
finite temperature. The choice of the model was made basically from the requirement of the simplicity of analysis. 
Certainly, the generalization to the case of long range gauge models would be of significant interest. 

Using the Hubbard-Stratonovich method of auxiliary fields, we sought for tachyonic excitations in the normal state 
of the NJL2+1 model in a magnetic field at finite temperature. We would like to note that the consideration of the 
theory at finite temperature is a necessary feature of our analysis. As discussed in Sec. II, although the symmetric 
state of the effective potential is unstable, its quadratic form of fluctuations that follows from ((7J is positive definite, 
hence tachyonic excitations are absent. The situation changes at finite temperature, where the effective potential (TTU)) 
has the instability typical for a second order phase transition. Utilizing the Schwinger-Keldysh real-time formalism, 
the dispersion relations for the composite fields ipip an d ipil 5 ^ were calculated in the LLL approximation in Sec. Ill, 
and for temperature less than a critical one a tachyonic excitation in the normal ground state was found. Thus, 
although there is no resonance state in the quantum mechanical one particle problem, the corresponding quantum 
field-theoretic problem in a magnetic field does have a tachyonic excitation in the normal state for temperature less a 
critical one. The contribution of higher Landau levels into dispersion relations for the composite fields a and 7r was 
taken into account in Sec. IV. 
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Appendix A: Closed form for the correlators U cr ' 7T 

To perform the summation over the Landau levels in Eq. (|55|) we assume that T < s/\eB\. Then wc can set tfj I — 1 
for all m > 1 while keeping to = tanh(a/2T), and expression (|55l) takes the form 



n^m k ) - ftanh -2|e5Hl±l)-» 2 (l + l) 

II (O, k) - wp (tanh ^ 1 j ^ _ + _ _ 



1 E m + En ( , . _ ±s nm (y)a 2 - 2\eB\r nm {y) 

ttZ 2 (E m + E n ) 2 - W { Snm[V) E n E m 

n,rn— 



Thus, in the considered approximation, the temperature dependence is described by the terms in the first line of the 
above equation. To calculate the first sum over the Landau levels in Eq. (|A1|) we use the representation 1/a = J °° dte~ at 
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valid for Re a > 0, and take into account that son(y) = Sn.o(y) = 2/™ e v /2n\. The evaluation of the second sum in 
Eq. ([Al[) is more involved. First, we use the chain of transformations 



(E n + E m ) 2 - VL 2 V ' E n E, { 



du {uj(u) + id), 1) 



-co 

oo 

1 f dhdh 2^L- tl E*-t 2 E 2 fh + t 2 + 2t 1 t 2 Q 2 

/ p *l+'2 " ••ml 1 

V*J Vh + h V z(h + t 2 ) 2 ' 



valid for Q 2 < E%. Then the sum 



n,m—0 

is evaluated using the integral representation ([52ft and the summation formula, 

oo 

^2 L«(z)x n = (1 - x)-( Q+1 ) exp 



n=0 



X- 1 



X < 1. 



Finally, the space integral over r in Eq. (l52l gives 

S(h,t 2 ) = icoth(|eB|(^ + i 2 ))exp (~a 2 (h +t 2 ) - 

Similarly, for another sum we obtain (r n o(y) = r 0m (y) = 0), 

e -ff 2 (t 1 +t 2 ) 



2y sinh |eB|ti sinh |eS|t2 
sinh|e£|(fi + t 2 ) 



R(h,t 2 )= r nm (y)e~ tlE -- 



'2-Ef, 



4sinrT \eB\(ti + t 2 ) 



exp 



2y sinh |eB|ti sinh |e£?|i 2 



sinh|eB|(ti +t 2 ) 



x 1- 



2y sinh |e.B|ii sinh |e£?|£ 2 



sinhleBlOi +i 2 ) 



Thus, we get the following representation for the correlators: 



IF'Vfi.fc) = -^ftanh— -1 
V ' ' Til 2 V 2T 



1 lTl 

4ct 2 - ft 2 2 



- / tft 1 dfc 2 e< tl+ta > ni, + 2 * n < tl -* a ) 



-2|eB|(t 1 +t 2 ) 1 ± 1 



n 2 (e ye 



-2|aB|(ti+ta) 1 =P 1 



1 /" d*idt 2 £^!lila Ai +i 2 + 2t!f 2 fi 2 =f2ct 2 (<i +f 2 )' 



2(ii+t 2 ) 5 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



1)- 



S(h,t 2 )+2\eB\R(t 1 ,t 2 )), (A7) 



which is convenient for expansions in k 2 and f2 2 . It is also very useful for obtaining the zero field limit, for that we 
get 



1 



2tt 3 / 2 



dpe~ s 



°- + p({n 2 -tf)x(l-x)Tv 2 ) 



(A8) 



l/A 2 



where an ultraviolet cutoff A is introduced at the lower limit of integral. 

It is obvious that the contribution of the first term in square brackets in Eq. (|A7j) is given by 



n?(o,A;) = o, n^(n,ifc) 



2 cr (tanh ^, - l) _ y 
ri 2 4ct 2 - Vt 2 6 ■ 



(A9) 
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The contribution of other terms can be expanded in y and f2 2 /|eP|, and keeping only the first order terms we get 



where 



and 



1-tanh^) +_^n-(fi,A:), 



ir 3 / 2 l 



n+(fi,A0, 



n 2 i 2 



-Q ± +yP ± , 



(A10) 
(All) 

(A12) 

(A13) 
(A14) 
(A15) 

The integral in the expression for V ± is divergent and we regularized it by introducing a lower limit cutoff e = 1/A 2 ^ 2 . 
Finally, we get 



dp 



(1 ± 2m 2 p) cothp 



2p 



sinh p 



2 -2;2 



- / dp Vpe- 



1 ± 



P± = 



VP 



1 ± 2m 2 p 
4p" 



2m 2 p 



coth p - 



coth p + 



1 



2p 



3 sinh p 



sinh p 



(1 — pcothp). 



p-> 
p- 



3 /tt 
V 2 

1 /¥ 
2V2 
^ 1 



, 1 m 
C (2' 1 + T 



to 2 \ 1 =F 1 1 ± 1 



2^ 



, 3 , m 2 



'5 m 2 \ 1T1 1±1' 
m C — , 1 H t= 



2 V2m' 

HF ( 2 71 m 2 \ m 4 /3 m 2 \ / - 



(A16) 
(A17) 
(A18) 
(A19) 



where C( s i w ) i s the generalized zeta function. 
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